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HIGHLIGHTS 


•  An  electro-thermal  model  capturing  battery  SOC,  voltage,  skin  and  core  temperature. 

•  A  convenient  parameterization  method  identifying  the  two  sub-models  separately. 

•  Validation  of  the  identified  model  with  electrochemical  impedance  spectroscopy. 

•  Experimental  validation  with  cycles  covering  20-100%  SOC,  5-38  °C,  and  max  C-rate  20C. 
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Combining  several  existing  lumped-parameter  models,  this  paper  presents  an  electro-thermal  model  for 
cylindrical  batteries.  The  model  consists  of  two  sub-models,  an  equivalent-circuit  electrical  model  and  a 
two-state  thermal  model  which  are  coupled  through  heat  generation  and  temperature  dependence  of 
the  electrical  parameters.  The  computationally  efficient  5-state  model  captures  the  state  of  charge  (SOC), 
terminal  voltage,  surface  temperature  and  the  often  neglected  core  temperature  of  a  battery  for  wide 
range  of  operating  conditions.  The  proposed  parameterization  scheme  allows  separate  identification  of 
the  electric  and  thermal  sub-models,  greatly  reducing  the  complexity  of  the  parameterization  process. 
The  methodology  is  applied  to  a  LiFeP04/graphite  battery.  Comparison  with  the  electrochemical 
impedance  spectroscopy  data  clarifies  the  frequency  range  of  the  model  fidelity.  The  model  is  further 
validated  with  two  drive-cycle  tests,  covering  SOC  range  25%— 100%,  temperature  5  °C— 38  °C,  and 
maximum  C-rate  of  22C. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Vehicle  electrification,  including  application  of  hybrid  electric 
vehicles  (HEV),  plug-in  electric  vehicles  (PHEV)  and  battery  electric 
vehicles  (BEV),  has  been  considered  as  a  promising  solution  to  the 
upcoming  energy  and  environmental  issues.  Among  all  types  of 
batteries,  lithium  ion  batteries  are  nowadays  widely  used  for 
automotive  applications  due  to  their  advantages  in  energy/power 
density,  charge  efficiency  among  others.  Correct  estimation  of  the 
electrical  and  thermal  dynamics  of  batteries  is  critical  for  safe  and 
efficient  operation.  Model  based  observers  have  been  widely 
considered  for  such  purpose,  where  electrical  models  are  used  to 
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estimate  the  state  of  charge  (SOC)  and  voltage  [1-3],  and  thermal 
models  are  applied  for  temperature  monitoring  [4-7], 

The  criteria  for  judging  the  quality  of  a  control-oriented  model 
include  model  fidelity,  computational  efficiency  and  ease  of 
parameterization.  Existing  electrical  models  vary  in  complexity  and 
fidelity.  In  some  applications,  where  the  applied  current  is  small 
and  simple,  a  coulomb  counting-based  OCV-R  model  is  sufficient 
[8],  Such  model  consists  of  an  SOC-dependent  open  circuit  voltage 
(OCV)  and  an  ohmic  resistance.  More  complicated  models,  such  as 
partial  differential  equation  (PDE)-based  electrochemical  models 
[9-11],  have  also  been  used  to  capture  the  electrochemical  pro¬ 
cesses  during  battery  operation,  giving  more  accurate  voltage 
estimation  over  a  wide  range  of  current  inputs.  However,  these 
models  are  computationally  intensive  and  difficult  to  parameterize 
[12],  To  balance  complexity  and  fidelity,  equivalent  circuit  models 
[1,13,14]  are  widely  adopted  in  battery  management  systems  (BMS), 
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where  a  series  of  resistor-capacitor  (R— C)  pairs  are  used  to  capture 
the  voltage  dynamics. 

For  temperature  estimation  of  the  cylindrical  battery,  existing 
PDE-based  thermal  models  can  predict  the  detailed  temperature 
distribution  throughout  the  battery  [10,15,16],  However,  like  the 
electrochemical  models,  these  models  are  not  easy  to  compute  in 
real  time  and  parameterize.  Therefore,  single-state  models,  where 
the  thermal  dynamics  are  represented  by  a  bulk  temperature 
[17,18],  have  been  widely  used  due  to  their  computational  effi¬ 
ciency.  Nevertheless,  such  approximation  might  lead  to  over¬ 
simplification  since  the  temperature  in  the  battery  core  can  be 
much  higher  than  in  the  surface  [19],  Two-state  thermal  models 
capturing  both  the  surface  and  the  core  temperatures  have  also 
been  studied  in  Refs.  [6,19—21],  Such  models  are  capable  of  esti¬ 
mating  the  core  temperature,  and  are  efficient  for  control  practice 
due  to  their  limited  number  of  states. 

Electro-thermal  models  constructed  by  combining  the  electrical 
and  thermal  models  introduced  above  have  also  been  studied 
broadly.  Some  of  these  works  are  conducted  with  pure  PDE-based 
approaches,  such  as  [15,22],  where  an  electrochemical  model  is 
coupled  with  a  heat-equation  based  thermal  model.  Some  apply 
lumped  models  to  capture  both  the  electrical  and  thermal  dy¬ 
namics,  such  as  [23]  (equivalent-circuit  electrical  model  +  single¬ 
state  thermal  model),  and  [4]  (equivalent-circuit  +  lumped  thermal 
model  along  the  longitudinal  direction).  Others  apply  mixed  ap¬ 
proaches.  For  example,  in  Refs.  [10,24],  the  electrical  dynamics  are 
modeled  by  an  electrochemical  model,  and  a  single-state  thermal 
model  is  used  to  capture  the  bulk  temperature. 

In  this  paper,  a  new  control-oriented  electro-thermal  model  is 
formulated  by  integrating  an  equivalent-circuit  electrical  model 
with  a  two-state  thermal  model.  The  thermal  model  captures  the 
surface  and  core  temperatures  as  the  two  states.  Compared  with  the 
single-state  model,  this  model  provides  higher  fidelity  with  mini¬ 
mal  increase  in  complexity,  and  it  is  more  computationally  efficient 
than  the  PDE-based  models.  Parameters  of  the  electrical  model  are 
dependent  on  temperature,  SOC  and  current  direction.  As  a  result, 
model  parameterization  can  be  computationally  expensive  and 
time  consuming,  requiring  large  dataset  and  advanced  optimization 
methods,  such  as  the  genetic  algorithm  [1[.  The  parameterization 
strategy  developed  in  this  paper  allows  separate  identification  of 
electrical  and  thermal  parameters  in  a  straightforward  way  by 
isolating  the  dynamics  of  each  sub-model.  The  electrical  parameters 
are  first  identified  under  isothermal  and  SOC-invariant  conditions 
for  different  temperatures,  SOCs  and  current  directions,  similar  to 
the  method  in  Ref.  [25],  The  thermal  model  is  then  parameterized 
with  a  drive-cycle  test  using  the  heat  generation  calculated  based 
on  the  current  and  voltage  measurement  and  the  open  circuit 
voltage  in  the  electrical  model.  The  designed  method  has  been 
applied  to  an  A123  LiFePC^/graphite  battery.  As  an  extension  of  [26], 
in  this  paper,  frequency  response  of  the  modeled  voltage  dynamics 
is  studied  by  comparing  the  model  impedance  with  the  electro¬ 
chemical  impedance  spectroscopy  (EIS)  measurement.  The  identi¬ 
fied  electro-thermal  model  is  also  validated  by  two  drive-cycle  tests. 
Good  model  fidelity  is  achieved  under  high  current,  shown  by  a 
drive-cycle  test  with  maximum  current  rate  of  22C  and  temperature 
variation  between  25  °C  and  38  °C. 

2.  Coupled  electro-thermal  battery  model 

In  this  section,  a  coupled  electro-thermal  model  is  formulated 
for  cylindrical  batteries.  The  terminal  voltage  is  captured  by  an 
equivalent  circuit  model,  and  a  two-state  thermal  model  is  adopted 
to  estimate  the  core  and  surface  temperatures.  The  parameters  of 
the  electrical  model  depend  on  temperature,  SOC  and  current  di¬ 
rection,  while  the  thermal  parameters  are  treated  as  constant. 


2.1.  Electrical  model 

The  schematic  of  an  equivalent  circuit  model  is  shown  in  Fig.  1. 
Under  a  current  /  (positive  for  discharge),  the  terminal  voltage  W  is 
modeled  in  three  parts  as 


VT  =  Vocv-IRs-±VRCi.  (1) 

i  =  1 

The  first  part  is  the  open  circuit  voltage,  Vocv,  which  is  usually  a 
nonlinear  function  of  SOC.  The  SOC  is  calculated  by  coulomb 
counting  as 


dSOC  _  % 

dt  “  Cbat’ 


(2) 


where  Cbat  is  the  battery  capacity.  The  second  part  is  the  voltage 
drop  over  a  series  resistor  Rs,  accounting  for  the  ohmic  resistance. 
The  last  part  is  the  voltage  drop  across  a  series  of  parallel  R— C 
circuits,  VRCii,  which  are  used  to  approximate  the  voltage  dynamics 
during  transients.  The  dynamics  of  each  R-C  pair,  Vrc.i,  are 
described  as 


dVRC,i 

dt 


(3) 


where  Rj  and  Q  are  the  equivalent  resistance  and  capacitance. 
Typically,  the  values  of  R,  and  Q  vary  with  SOC,  current  direction, 
and  temperature  which  is  captured  by  the  thermal  model  to  be 
introduced  in  Section  2.2.  It  will  be  shown  in  Section  3.1  that  in¬ 
clusion  of  more  R-C  pairs  improves  the  model  fidelity  at  the  cost  of 
complexity. 


2.2.  Two-state  thermal  model 


By  assuming  longitudinal  homogeneity,  a  two-state  model 
[6,20],  as  shown  in  Fig.  2,  is  used  to  capture  the  lumped  thermal 
dynamics  of  a  cylindrical  battery.  Governing  equations  of  Tc  and  Ts 
are  [6,20], 

Cc§  =  0-  +  ^ 

(4) 

Cs§  =  2jp- 

The  heat  generation  Q  in  Eq.  (4)  is  a  byproduct  of  the  chemical 
reactions  taking  place  in  the  electrode  assembly  during  battery 
operation.  The  value  of  Q  is  calculated  based  on  the  electrical 
model,  as 

Q  =  I(V ocv  -  Vr),  (5) 


Rl  Rn 
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Coolant  Convection 


Fig.  2.  Schematic  of  the  two-state  thermal  model  (adapted  from  Ref.  [20]). 


which  includes  the  joule  heating  and  energy  dissipated  in  the 
electrode  over-potentials  [18],  For  more  accurate  modeling  of  heat 
generation,  reversible  heat,  such  as  entropic  heat  can  be  included 
[10,19],  A  heat  conduction  resistance,  Rc,  is  used  to  model  the  heat 
exchange  between  the  core  and  the  surface.  Convective  cooling 
through  the  battery  surface  is  modeled  by  a  convection  resistance, 
Ru,  between  Ts  and  Tf.  The  value  of  Ru  depends  on  the  pack  geom¬ 
etry,  type  of  coolant  and  coolant  flow  rate.  The  rates-of-change  of  Tc 
and  Ts  are  determined  by  the  heat  capacity  of  the  core,  Cc,  and  that 
of  the  casing,  Cs. 


2.3.  Electro-thermal  coupling 

The  electrical  and  the  thermal  sub-models  interact  with  a  two- 
way  coupling,  as  shown  in  Fig.  3.  First,  the  electrical  model  calcu¬ 
lates  the  SOC  and  voltage  of  the  battery  based  on  the  current  /  and 
the  electrical  parameters  Rs,  R,  and  Q.  The  difference  between  Vi¬ 
and  1/ocv,  along  with  /,  is  used  to  determine  the  heat  generation  Q. 
The  thermal  model  then  calculates  Tc  and  Ts  based  on  Q  and  the 
ambient  temperature  Tf.  The  core  temperature,  Tc,  which  represents 
the  temperature  of  the  lumped  electrode  assembly,  is  used  to 
determine  the  temperature-dependent  parameters  of  the  electrical 
model.  The  electro-thermal  coupling  complicates  the  parameteri¬ 
zation  of  the  model,  which  is  usually  conducted  by  fitting  the  pa¬ 
rameters  to  match  the  measured  outputs  (Vp  and  Ts)  from 
experiments.  If  the  temperature  is  changing  during  the  experi¬ 
ments,  the  temperature-dependent  electrical  parameters  will  also 
be  varying,  making  them  harder  to  identify.  Experiments  designed 
to  excite  the  voltage  dynamics  without  significantly  changing  the 
SOC  or  temperature  are  shown  in  Section  3.  Therefore,  the  pa¬ 
rameters  will  be  constant  (in  each  step  of  the  experiments),  and 
easy  to  identify. 


Fig.  3.  Coupling  between  the  electrical  and  the  thermal  models. 


Time  (min) 

Fig.  4.  Voltage  and  current  profile  during  the  discharging  pulse  test. 


3.  Model  parameterization 

In  this  section  the  methodology  for  parameterizing  the  electrical 
and  thermal  sub-models  is  described.  First,  the  electrical  parame¬ 
ters  are  identified  based  on  a  series  of  pulse-relaxation  tests.  Each 
test  is  performed  under  isothermal  (and  SOC-invariant)  conditions. 
In  this  way,  the  electrical  model  is  decoupled  from  the  thermal 
model  and  can  be  identified  independently.  The  tests  are  repeated 
under  different  SOCs,  temperatures  and  current  directions.  The 
thermal  model  is  then  parameterized  with  the  heat  generation 
calculated  based  on  the  measured  current,  voltage  and  the  identi¬ 
fied  open  circuit  voltage  in  the  electrical  model. 

3.1.  Parameterization  of  the  electrical  model 

Parameterization  of  the  electrical  model  involves  identification 
of  the  open  circuit  voltage  Vocv,  capacity  Cbat.  ohmic  resistance  Rs, 
and  resistance  and  capacitance  of  the  R-C  pairs. 

First,  Vocv  and  Cbat  are  determined  by  cycling  the  battery  at  a 
low  current,  e.g.  C/20.  The  battery  is  charged  with  the  constant 
current  (CC)  to  the  high  voltage  limit  recommended  by  the 
manufacturer  before  switching  to  constant  voltage  (CV)  charging 
until  the  current  drops  to  the  cut-off  limit,  e.g.  C/50.  The  battery  is 
then  discharged  to  the  low  voltage  limit  with  the  same  constant 
current  used  in  charge.  Since  the  voltage  drop  across  the  battery 
internal  resistance  is  negligible  under  the  low  current,  the  differ¬ 
ence  between  the  charging  and  the  discharging  voltage  is  usually 
small.  Therefore,  either  the  charging  or  discharging  voltage  profile 
can  be  used  as  Vocv-  However,  for  some  battery  chemistries,  e.g. 
LiFeP04,  the  difference  between  the  charging  and  discharging 
voltage  can  be  much  larger  due  to  hysteresis  [27,28],  In  such  case, 
the  average  of  the  charging  and  discharging  voltage  can  be  taken  as 
the  open  circuit  voltage,  and  the  hysteresis  in  voltage  may  need  to 
be  modeled.  The  100%  SOC  is  defined  at  the  end  of  the  CV  charge 
period,  and  0%  at  the  end  of  the  discharge.  The  capacity  Cbat  is 
determined  by  the  amount  of  charge  between  0%  and  100%  SOCs. 

The  parameters  of  the  equivalent  circuits,  including  Rs,  R;  and  Q, 
are  identified  based  on  the  pulse-relaxation  experiments.  Tests  are 
conducted  in  both  discharge  and  charge  current  directions,  and  the 
discharge  case  is  taken  as  an  example  for  illustration.  The  battery  is 
first  charged  to  100%  SOC  under  the  CC— CV  protocol  and  then 
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relaxed  for  2  h  to  reach  electrical  and  thermal  equilibrium.  It  is  then 
discharged  to  90%  SOC  with  a  current  of  1  C  (in  around  6  min), 
which  elevates  the  battery  temperature  only  minimally.  The  current 
is  then  cut  to  relax  the  battery  for  2  h.  The  voltage  evolution  during 
relaxation  is  recorded  for  identifying  the  R— C  parameters.  The 
procedures  are  repeated  until  the  voltage  reaches  the  low  voltage 
limit.  During  the  experiment,  the  battery  is  placed  inside  a  thermal 
chamber,  where  the  temperature  is  controlled  at  the  desired  set- 
point.  To  obtain  the  temperature  dependency  of  the  parameters,  the 
test  is  performed  at  multiple  temperatures.  The  voltage  and  current 
profiles  during  the  discharge  pulse  test  of  a  battery  are  plotted  in 
Fig.  4.  Voltage  relaxation  after  each  pulse,  as  shown  in  the  inset  of 
Fig.  4,  is  used  to  determine  the  parameters  of  the  equivalent  circuit 
at  each  SOC.  It  is  noted  that  since  both  the  temperature  and  SOC  are 
constant  during  the  relaxation,  the  parameters  can  be  easily  iden¬ 
tified  at  fixed  temperature.  The  ohmic  resistance  Rs  is  calculated  by 
the  initial  voltage  jump  AVS  based  on  the  Ohm’s  law, 


In  the  form  of  Eq.  (9),  we  have 
observation  :  Z  =  s2Ts(s) 

regressors  :  4>  =  [Q(s),  Tf(s)-Ts(s),  sTs(s)]t  (11) 

parameters  :  ft  —  ![  a,  (H,  y]T, 

where  the  original  physical  parameters  are  now  lumped  into 


(12) 


The  least  squares  algorithm  is  applied  to  identify  8  [30]  based  on 
the  parametric  model, 


Rs 


AVS 


w  s  -  (*t*)~vz. 


(13) 


where  /  is  the  current  applied  during  the  pulse.  The  values  of  R{s 
and  Q’s  are  identified  by  using  the  relaxation  voltage  after  the 
initial  jump,  Vrex(tr),  shown  in  Fig.  4.  The  evolution  of  Vrex  can  be 
derived  by  solving  Eq.  (3),  as 

Vrex(tr)  =  IR'  (]  “  e~^j  (l  -  e  &) ,  (7) 

where  tp  is  the  pulse  duration  before  the  relaxation.  The  resistance 
Ri  and  the  capacitance  Q  of  each  R-C  pair  are  then  fitted  to  mini¬ 
mize  the  least  square  error  in  voltage, 

J  =  E  OWfc)  -  WataW)2,  (8) 


where  6  is  the  estimate  of  6.  The  parameter  identification  is  per¬ 
formed  in  the  time  domain,  and  the  matrices  Z  and  4>  consist  of  the 
signals  Z  and  $  at  different  time  instant  tk  during  the  experiment, 


r*T(t,n 

0T(fe) 

z  = 

,  $  = 

"W 

0T(tN) 

LmiSyJ 

m(tk)  =  \J\  +  0T(tk)0(tk), 


where  Vrex,data(k)  is  the  measured  relaxation  voltage  at  each  time 
instant  k.  The  number  of  R— C  pairs  (n)  can  be  determined  based  on 
the  fitting  results  to  be  discussed  later. 


3.2.  Parameterization  of  the  thermal  model 


The  parameters  of  the  thermal  model,  including  the  conduction 
resistance  Rc,  convection  resistance  Ru,  and  heat  capacities  of  the 
battery  core  and  surface  Cc  and  Cs,  are  identified  based  on  the 
method  developed  in  Refs.  [6,29],  Considering  the  difficulty  of 
measuring  the  core  temperature  Tc,  parameterization  is  designed 
by  only  using  the  current  I,  surface  temperature  Ts  and  coolant  flow 
temperature  If,  which  can  all  be  easily  measured. 

As  the  first  step,  a  parametric  model  in  the  Laplace  domain 
Refs.  [30], 

Z(s)  =  0T$(s)  (9) 


needs  to  be  derived.  In  Eq.  (9),  Z  is  referred  to  as  the  observation,  8  is 
the  parameter,  and  <I>  is  the  regressor  consisting  of  measured  sig¬ 
nals.  For  the  thermal  model  in  Eq.  (4),  the  parametric  model  can  be 
derived  by  performing  Laplace  transform  to  the  governing  equa¬ 
tions  and  eliminating  the  unmeasured  Tc  by  replacing  with 
measured  signals,  as 


+  (T'w  - 


Tc  +  Cs  1  N 

,CcQRc+QRu, 


)sTs(s). 


(10) 


where  z  and  0  are  the  time-domain  signals  of  Z  and  <I>,  m  is  the 
normalization  factor  used  to  enhance  the  robustness  of  the  algo¬ 
rithm,  and  N  is  the  number  of  time  instants.  In  Eq.  (14),  the  needed 
time-domain  signals  z  and  0  include  Q,  Ts,  Tf,  and  the  first  and 
second  time  derivatives  of  Ts  ( sTs  and  s2Ts  in  the  time  domain).  The 
temperature  Ts  and  Tf  are  measured,  and  Q,  is  calculated  based  on 
Eq.  (5).  As  for  the  time  derivatives  of  Ts,  in  order  to  avoid  differ¬ 
entiating  the  signals  which  is  easily  corrupted  by  noises,  a  second 
order  filter 


1  _  1 

A(s)  (b + Xj){S 


(15) 


needs  to  be  used.  The  filtered  parametric  model  will  become 


Zf 

Zf 


0T<f>  f 
1 

A(s) 

1 

W) 


(PT.)  , 

[Q  Tf- 


Jig 

A (S)  is 
Ts  STS]T 


[A®  «(»-'>■•)  iwT>] 


(16) 


The  time-domain  signals  ofZf  and  4>f,  Zf  and  <pf  will  be  used  in  Eq. 
(14)  in  place  of  z  and  0,  which  are  obtained  by  passing  the 
measured/calculated  signals  through  relevant  filters,  as  shown  in 
Fig.  5. 

The  95%  confidence  interval  of  each  identified  parameter  8j  can 
be  calculated  as  [31] 
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SOC 


Fig.  5.  Filtering  of  the  signals  used  to  identify  the  Thermal  Model. 


0j  ± 96  Jvartdj 


(17) 


The  estimated  variance  of  6j,  var(ffj),  is  the  diagonal  element  of 
the  covariance  matrix, 


cov(Tj  -m  ](VV)  V,  (18) 

where  a2  is  the  variance  of  Zf, 

„2  _  £„(*#„) -?V„))2  (|9) 

By  using  the  above  algorithm,  the  lumped  parameters  a,  0  and  y 
are  identified.  It  is  noted  that  only  three  of  the  four  original  pa¬ 
rameters  can  be  obtained  based  on  a,  /?,  and  y.  Therefore,  one  of  the 
parameters  needs  to  be  assumed  in  order  to  determine  the  other 
three.  It  is  found  that  Cs  can  be  easily  calculated  based  on  the 
specific  heat  capacity  and  dimensions  of  the  casing.  The  estimates 
of  Ru,  Rc  and  Cc  can  then  obtained  by 


_jCsRu  +  SC2R  u  +  1 

2r„cs 


(20) 


4.  Experiment  setup 

Experiments  have  been  conducted  on  a  2.3  Ah  A123  26650 
LiFeP04  battery  to  parameterize  and  validate  the  electro-thermal 
model.  A  Yokogawa  GS-610  Source-Measure  Unit  was  used  for  the 
tests  to  parameterize  the  electrical  model.  The  C/20  CC— CV  test  was 
performed  first  to  identify  Vocv  and  Cbat-  The  high  and  low  voltage 
limits  for  this  battery  are  3.6  V  and  2  V,  and  the  cutoff  current  of  CV 
is  set  at  50  mA.  The  unit  was  then  programmed  to  run  the  pulse- 
relaxation  tests  to  parameterize  the  R— C  pairs.  During  the  tests, 
the  battery  was  cycled  inside  a  Cincinnati  Sub-Zero  ZPHS16-3.5- 
SCT/AC  thermal  chamber,  where  the  ambient  temperature  was 
controlled.  The  test  is  conducted  between  5  °C  and  45  °C  at  every 
10  °C  interval  to  obtain  the  temperature  dependency  of  the 
parameters. 

A  drive-cycle,  the  scaled  urban  assault  cycle  (UAC)  [32],  which 
contains  sufficient  excitation  [6],  was  used  to  parameterize  the 
thermal  model.  The  test  was  conducted  with  a  Bitrode  FTV1-200/ 
50/2-60  cycler,  which  is  capable  of  providing  the  power  needed 
to  simulate  the  drive-cycle.  The  battery  was  placed  in  a  designed 
flow  chamber  with  a  fan  regulating  the  air  flow  rate.  The  fan  speed 
was  fixed  so  that  a  constant  convection  resistance  Ru  would  be 
identified.  T-type  thermocouples  were  installed  on  the  battery 
surface  to  measure  Ts,  and  inside  the  flow  chamber  to  measure  Tf. 
The  measured  Ts  and  Tf  were  used  for  model  identification.  In  order 
to  validate  the  identified  thermal  model,  a  thermocouple  was  also 
installed  inside  the  battery  to  measure  Tc  [6],  The  whole  experiment 
setup  is  shown  in  Fig.  6. 


Fig.  6.  Experiment  setup. 
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Fig.  8.  Rs  under  different  temperatures,  SOCs  and  current  directions. 


5.  Parameterization  results 

In  this  section,  the  parameters  of  the  LiFeP04/graphite  battery 
identified  based  on  the  methodology  discussed  in  Section  3  are 
demonstrated. 


SOC  over  the  tested  range  (10%-90%  SOC),  moderate  dependence 
on  current  direction,  and  strong  dependence  on  temperature.  The 
value  of  Rs  decreases  as  the  temperature  rises,  dropping  by  more 
than  50%  from  19.1  mQ  at  5  °C  to  8.2  mO  at  45  °C  under  discharge. 
This  trend  can  be  explained  by  the  more  active  reaction  kinetics  at 
higher  temperature.  An  Arrhenius-like  function, 


5.1.  Identified  parameters 

The  open  circuit  voltage  Vocv  and  capacity  Cbat  are  first 
measured  based  on  the  C/20  cycling  test.  The  measured  terminal 
voltage  is  shown  in  Fig.  7.  As  mentioned,  due  to  the  hysteresis  of  the 
LiFePC>4  chemistry,  a  gap  is  observed  between  the  charge  and 
discharge  voltage  curves.  The  average  of  the  two  is  taken  as  Vocv 
and  modeling  of  the  hysteresis  remains  to  be  addressed.  The  actual 
battery  capacity  Cbat  was  measured  at  2.1  Ah. 

The  series  resistance  Rs  identified  from  the  pulse-relaxation 
tests  under  different  SOCs,  temperatures  and  current  directions  is 
shown  in  Fig.  8.  It  can  be  seen  that  Rs  has  a  minimal  dependence  on 


(21) 


is  used  to  fit  the  temperature  dependency  of  Rs,  where  /  =  c  for 
charge,  and  1  =  d  for  discharge. 

To  parameterize  the  rest  of  the  electrical  model,  the  number  of 
R— C  pairs  needs  to  be  determined  first.  Equivalent  circuits  with 
single,  double  and  triple  R-C  pairs  have  been  fitted  to  the  relaxa¬ 
tion  voltage  in  the  pulse-relaxation  test.  Comparison  of  the  fitting  is 
shown  in  Fig.  9.  It  can  be  seen  that  the  single  R-C  model  yields  large 
error  especially  during  the  first  500  s  of  relaxation,  whereas  the 
error  of  the  double  R— C  model  is  much  smaller  over  the  entire 
period.  The  improvement  from  the  double  R-C  to  the  triple  R— C 
model  is  much  smaller.  The  double  R-C  model  is  hence  chosen  for 
its  advantage  in  computational  efficiency.  The  identified  parame¬ 
ters  of  the  double  R— C  pairs  are  shown  in  Figs.  10—13. 

For  Ri  and  R2,  it  can  be  seen  that,  similar  to  Rs,  their  values 
decrease  as  the  temperature  increases  due  to  the  more  active  ki¬ 
netics  at  elevated  temperature.  Nevertheless,  compared  with  Rs,  Ri 
and  R2  exhibit  more  prominent  dependence  on  SOC.  Interestingly, 
such  dependence  is  opposite  for  charge  and  discharge.  The  charge 
resistance  increases  significantly  in  the  high  SOC  range  (80%- 
100%),  whereas  the  discharge  resistance  is  higher  at  SOCs  near  0%. 
The  observed  SOC  dependency  is  in  accordance  with  those  in 
Refs.  [25],  where  the  parameters  identified  at  SOCs  closer  to  0%  and 
100%  are  available.  A  second  order  polynomial  is  used  to  fit  the  SOC- 
dependency  of  R]  and  R2,  and  an  Arrhenius-like  function  for  tem¬ 
perature  dependency, 


Ri:i  =  (%i,i  +  Qi.i.iSOC  +  a2,uSOC2)  -exp^^J : 

f  1,  for!?!  ^  ( c,  for  charge  (22) 

\  2,  for  R2  \  d,  for  discharge  ’ 


Fig.  9.  Fitting  of  the  relaxation  voltage. 


Fig.  12.  C, 


■  different  temperatures,  SOCs  and  current  directions. 
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where  Ro,i,i,  am,  and  024,1  are  the  fitted  parameters  for  the  SOC 
dependency,  and  TrefR,i,i  and  TShjftR,i,i  are  for  the  temperature  de¬ 
pendency.  The  identified  capacitances  C\  and  C2,  plotted  in  Figs.  12 
and  13,  also  show  notable  dependence  on  temperature  and  SOC.  As 
opposed  to  the  resistance,  both  capacitance  values  increase  with 
temperature.  For  the  SOC  dependence,  Ci  shows  a  nearly  mono- 
tonically  increasing  trend  from  low  to  high  SOC  under  discharge, 
while  its  values  under  charge  are  relatively  constant  in  the  low  and 
middle  SOC  range  but  see  a  sharp  decrease  in  the  high  SOC  end.  The 
relationship  between  C2  and  SOC  is  more  complicated,  as  double 
peaks  are  observed.  The  values  of  Ci  and  C2  are  each  fitted  by  a 
quadratic  function  of  SOC  affine  with  temperature, 

0,1  =  C044  +  cul(SOC)  +  c2il(SOC)2  +  (c3iU  +  c4il  (SOC) 

+  c5,i,l(SOC)2)  Tm, 
l  _  /  1  >  for  i  _  /  c’  charge 
1 2,  for  R2  ’  \d,  discharge ' 

(23) 
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Fig.  14.  Measured  current,  battery  surface  and  flow  temperatures  under  UAC. 


The  thermal  model  is  parameterized  by  using  data  from  the  UAC 
drive-cycle  test.  The  measured  current  I,  coolant  flow  temperature 
Tf,  and  battery  surface  temperature  Ts  are  shown  in  Fig.  14.  The 
current  I  and  measured  V  are  used  to  calculate  the  heat  generation 
along  with  the  parameterized  open  circuit  voltage  based  on  Eq.  (5). 
The  pre-calculated  Cs  and  identified  Cc,  Rc  and  Ru,  along  with  the 
95%  confidence  interval,  are  listed  in  Table  1.  The  identified  thermal 
parameters  are  very  close  to  those  in  Refs.  [6],  where  the  heat 
generation  is  calculated  by  using  a  single  lumped  internal  resis¬ 
tance  Re,  as  Q  =  f2Re.  This  lumped  Re  includes  the  heat  dissipation 
through  the  series  resistance  Rs  and  the  R-C  pairs  of  the  equivalent 
circuit  model  in  this  paper.  As  pointed  out  in  Refs.  [6],  the  value  of 
Re  can  be  cycle-dependent,  or  even  non-constant  within  a  single 
cycle  due  to  the  variability  of  battery  resistance  on  temperature, 
SOC  and  other  factors.  Consequently,  if  a  single  Re  is  used,  its  value 
needs  to  be  obtained  by  online  identification  as  in  Ref.  [6],  The 
identified  Re  under  the  UAC  is  11.4  mO,  which  is  in-between  the 
series  resistance  Rs  and  the  total  resistance  Rs  +  Ri  +  R2  identified  in 
this  paper. 

6.  Validation  of  the  identified  model 

This  section  is  devoted  to  validating  the  identified  electro¬ 
thermal  model.  The  model  dynamics  are  first  investigated  by 
comparing  the  model  impedance  with  EIS  measurements.  The  er¬ 
rors  in  voltage  and  temperature  estimation  under  two  drive-cycle 
tests  will  also  be  evaluated. 


6.1.  Comparison  of  the  modeled  electrical  impedance  with  EIS 

The  battery  impedance,  Z,  is  the  relationship  between  the  cur¬ 
rent  I  and  the  terminal  voltage  Vt  as  a  function  of  the  current 
frequency, 


Table  1 

Identified  parameters  of  the  thermal  model. 


Value 


Q0K-1)  4.5 

Cc  CJK-1)  62.7 

RC(KW')  1.94 

Ru  (KW-1)  3.19 


95%  Confidence  interval 


59.1-66.4 

1.86-2.01 

3.19-3.20 
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Vj(j<o)  =  Z(jb))I(ju),  (24) 

where  j  is  the  imaginary  unit,  and  w  is  the  angular  frequency.  For 
the  electrical  model  in  Eqs.  (1)— (3),  the  impedance  at  To  and  SOCo, 
can  be  obtained  based  on  the  transfer  function  from  the  current  to 
the  voltage,  with  the  model  linearized  at  To  and  SOCo.  The  transfer 
function,  Z(s),  is  derived  by  taking  Laplace  transformation  of  the 
linearized  model,  as 

Vj(s)  =  Z(s)I(s) 

Z(s )  =  - — - - y  ,  (25) 

CbatS  ci(s  +  ro)  c2(s  +  r2g) 

where  a  is  the  slope  of  Vocv,  and  all  parameters  are  evaluated  at  To 
and  SOCo.  By  taking  s  =  jio,  the  impedance  Z(/w)  is  obtained.  The 
Nyquist  plot  of  the  modeled  Z,  where  the  real  part  RE(Z )  is  plotted 
against  the  imaginary  part  /M(Z),  at  25  °C  and  80%  SOC,  is  shown  in 
Fig.  15  as  an  example.  It  can  be  seen  that  the  modeled  impedance 
consists  of  two  arcs  in  the  high  and  middle  frequency  ranges,  cor¬ 
responding  to  the  dynamics  of  the  two  R— C  pairs.  The  intersection 
with  the  real  axis  represents  the  series  resistance  Rs. 

The  actual  battery  impedance  can  be  measured  by  EIS,  [33,34], 
where  the  battery  is  excited  by  AC  currents  with  different  fre¬ 
quencies  and  the  response  of  the  voltage  is  recorded  to  calculate 
the  impedance.  EIS  has  been  performed  under  the  open  circuit 
condition,  with  an  Autolab  PGSTAT302N  potentiostat/galvanostat 
under  50  mA  AC  currents,  sweeping  frequency  from  0.1  mHz  to 
5  kHz.  The  Nyquist  plot  of  the  measured  impedance  is  shown  in 
Fig.  15,  and  the  magnitude  and  phase  of  the  impedance  at  different 
frequencies  is  plotted  in  Fig.  16  to  compare  with  the  model.  As  seen 
from  Fig.  15,  the  measured  impedance  is  composed  of  three  parts,  a 
semicircle  in  the  frequency  range  higher  than  1  Hz,  corresponding 
to  the  dynamics  of  the  passive  films  at  the  solid/electrolyte  inter¬ 
face  (SEI),  an  arc  in  the  frequency  range  2  mHz— 1  Hz,  representing 
the  dynamics  of  the  charge-transfer  resistance  and  the  double  layer 
capacitance,  and  a  straight  line  below  2  mHz,  related  to  the  dy¬ 
namics  of  lithium  diffusion  [33,35]. 

The  modeled  impedance  is  more  accurate  between  0.01  Hz  and 
100  Hz  than  in  the  higher  and  lower  frequency  ranges,  according  to 
Fig.  16.  The  obtained  accuracy  in  impedance  is  comparable  to  that  in 
Refs.  [36]  and  [37],  where  similar  electrical  models  are  analyzed.  EIS 
has  also  been  conducted  at  other  SOCs,  such  as  60%  and  95%, 


RE(Z)  (n) 


Fig.  15.  Nyquist  plot  of  the  impedance:  model  vs.  EIS  measurement. 
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Fig.  16.  Magnitude  and  phase  of  the  impedance  at  different  frequencies:  model  vs.  EIS 
measurement. 


revealing  similar  behavior.  At  the  high  frequency  end,  as  shown  in 
Fig.  15,  the  small  semi-circle  representing  the  dynamics  of  the  SEI 
film  is  not  modeled,  as  the  right  end  of  the  semi-circle  is  captured 
as  the  series  resistance  Rs.  The  reason  may  be  attributed  to  lack  of 
fast  voltage  dynamics  excited  in  the  pulse-relaxation  tests.  In  the 
low  frequency  range,  the  impedance  is  underestimated  due  to  two 
possible  reasons.  As  indicated  by  Eq.  (25),  the  modeled  impedance 
at  low  frequencies  is  related  to  the  slope  of  Vocv, «.  as  well  as  R2  and 
C2.  For  one  thing,  a  may  be  inaccurate  because  of  the  flatness  of 
Vocv  (as  shown  in  Fig.  7)  and  exclusion  of  hysteresis  in  1/ocv 
modeling.  For  another,  the  pulse-relaxation  voltage  used  to  iden¬ 
tified  the  R-C  pairs  is  not  very  sensitive  to  the  slow  dynamics  ( R2 
and  C2),  as  compared  to  the  fast  dynamics  (Ri  and  Ci),  which  might 
compromise  the  accuracy  of  the  identified  R2  and  C2.  The  full  bat¬ 
tery  dynamics  can  be  modeled  more  accurately  based  on  the  EIS 
data  by  using  complex  impedance  components,  such  as  constant 
phase  elements  (CPE)  [33,34],  However,  the  modeling  and  param¬ 
eterization  procedure  will  be  much  more  complicated,  and  the 
implementation  of  model  in  time  domain  is  usually  achieved  by 
approximation  with  a  significant  number  of  R— C  pairs  [33],  The 


Fig.  17.  Current  profile  and  SOC  variation  of  the  charge-depleting  cycle  (top:  current  in 
time;  middle:SOC;  bottom:  FFT  of  current). 
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associated  computational  load  may  cause  issues  for  onboard  ap¬ 
plications,  especially  when  the  number  of  cells  are  large. 

6.2.  Validation  by  drive-cycles 

The  model  has  also  been  validated  by  data  from  two  drive-cycle 
tests.  The  current,  SOC,  and  temperature  variation  of  the  validation 
cycles  are  plotted  in  Figs.  17-20.  The  first  validation  cycle  is  a 
charge-depleting  (CD)  cycle  with  maximum  current  around  10  C 
and  average  at  4  C.  The  battery  SOC  decreases  from  100%  to  25%, 
and  the  core  temperature  increases  from  5  °C  to  23  °C  in  approxi¬ 
mately  10  min.  The  second  validation  cycle,  shown  in  Figs.  19  and 
20,  is  a  more  drastic  and  charge-sustaining  (CS)  cycle,  where  the 
maximum  current  is  more  than  20  C,  and  the  average  around  0.03  C. 
The  resulting  SOC  varies  between  30%  and  53%,  while  the  tem¬ 
perature  fluctuates  between  25  °C  and  38  °C.  The  two  cycles 
combined  cover  an  SOC  range  from  25%  to  100%,  and  a  temperature 
range  between  5  °C  and  38  °C,  which  are  roughly  the  ranges  that 
the  model  parameters  are  identified  in  this  paper.  Fourier  trans¬ 
formation  has  also  been  performed  to  demonstrate  the  frequency 
spectrum  of  the  two  cycles,  as  shown  in  Figs.  17  and  19.  It  can  be 
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seen  that  the  CD  cycle  is  slower  and  contains  a  peak  at  0  Hz,  which 
is  the  DC  component  corresponding  to  the  mean  charge-depleting 
rate.  The  dynamics  of  the  CS  cycle  are  faster  and  have  more  evenly 
distributed  frequency  content. 

The  voltage  and  temperatures  estimated  by  the  parameterized 
model  are  plotted  in  Figs.  18  and  20  to  compare  with  the  mea¬ 
surements.  The  root  mean  square  errors  (RMSE)  of  the  estimation 
by  using  parameters  directly  identified  and  based  on  fitted  func¬ 
tions  Eqs.  (22)  and  (23)  are  listed  in  Table  2.  It  can  be  seen  that  the 
identified  model  shows  good  fidelity  under  the  high-power  charge- 
sustaining  cycle,  as  the  voltage  RMSE  is  below  20  mV,  and  the 
temperature  RMSE  below  1°C.  Larger  voltage  RMSE  (45  mV)  are 
observed  in  the  charge-depleting  cycle.  The  model  performance  is 
compared  to  that  of  similar  electrical  models  under  drive-cycle 
tests  in  Refs.  [1,14],  In  Ref.  [1],  an  equivalent  circuit  electrical 
model  with  2  RC  pairs  and  1  hysteresis  state  is  parameterized  and 
validated  with  a  specific  drive-cycle  test.  In  Ref.  [14],  electrical 
models  with  1  -3  RC  pairs  and  1  hysteresis  state  are  validated  under 
the  Dynamic  Stress  Test  (DST)  and  Federal  Urban  Dynamic  Schedule 
(FUDS)  respectively.  The  specifications  of  the  models  and  drive- 
cycles  and  the  estimation  errors  are  shown  in  Table  3.  The 
voltage  RMSEs  of  models  in  Ref.  [14]  are  as  small  as  8  mV,  partially 
due  to  the  mild  operating  condition  (small  current  rate  and 
invariant  temperature)  besides  the  well  identified  model.  The 
estimation  error  of  the  electro-thermal  model  in  this  paper  under 
the  CS  cycle  is  close  to  that  of  the  model  in  Ref.  [1  ].  As  for  the  CD 
cycle,  there  are  several  possible  reasons  for  the  larger  voltage 
estimation  error.  First,  hysteresis  of  the  LiFeP04  chemistry  has  not 
been  included  in  the  modeling.  The  nonlinear  hysteresis  dynamics 
in  voltage  are  more  prominent  under  unidirectional  current  than 
under  charge-balance  current  profile.  Second,  it  has  been  revealed 
by  the  EIS  measurement  that  the  battery  voltage  dynamics  are  not 
very  well  captured  in  the  low  frequency  range  (below  0.01  Hz),  and 
the  mismatch  becomes  larger  as  the  frequency  approaches  0  Hz 
(Figs.  15  and  16).  The  current  profile  of  the  CD  cycle,  whose  fre¬ 
quency  distribution  has  a  high  peak  at  0  Hz  corresponding  to  an 


Fig.  19.  Current  profile  and  SOC  variation  of  charge-sustaining  cycle  (top:  current  in 
time;  middle:SOC;  bottom:  FFT  of  current). 


Table  2 

Estimation  errors  when  using  parameters  directly  identified  and  parameters  based 


RMSE  VT  (mV)  Ts  (°C)  Tc  (°C) 

CD  cycle  (directly  identified)  45.4  0.20  0.56 

CD  cycle  (fitted  function)  48.6  0.25  0.62 

CS  cycle  (directly  identified)  20.3  0.65  0.92 

CS  cycle  (fitted  function)  19.2  0.67  0.96 
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7.  Conclusion 

In  this  paper,  an  electro-thermal  model  for  cylindrical  batteries 
is  proposed  based  on  an  equivalent-circuit  electrical  model  and  a 
two-state  thermal  model.  The  electrical  and  the  thermal  parts  are 
coupled  through  heat  generation  and  temperature  dependency  of 
the  electrical  parameters.  The  integrated  5-state  model  is 
computationally  efficient  for  onboard  applications  and  provides 
estimates  of  the  critical  battery  states.  The  core  temperature  of 
cylindrical  batteries,  which  is  usually  neglected  by  existing 
control-oriented  models,  is  captured.  As  shown  in  Figs.  18  and  20, 
the  core  temperature  can  be  significantly  higher  than  the  surface 
temperature,  especially  under  high  current.  Accurate  estimation  of 
the  core  temperature  can  help  reduce  the  risk  of  overheating. 
Additionally,  since  the  parameters  of  the  electrical  model  are 
temperature  dependent,  capturing  the  temperature  close  to  the 
electrode  assembly  (core  temperature)  can  also  enhance  the  ac¬ 
curacy  of  voltage  estimation.  Such  advantage  is  desirable  espe¬ 
cially  under  high  current  application  where  the  temperature 
swing  is  large  and  fast. 

The  electrical  sub-model  is  parameterized  under  isothermal 
conditions.  Since  the  electrical  and  the  thermal  dynamics  are 
decoupled,  the  parameters  can  be  identified  separately.  Pulse- 
relaxation  tests  are  used  to  parameterize  the  electrical  model  so 
that  the  SOC  dependency  of  the  parameters  can  be  conveniently 
identified  under  SOC  invariant  conditions.  The  thermal  model  is 
then  parameterized  by  using  the  heat  generation  calculated  by  the 
measured  current,  voltage  and  the  electrical  model.  Such  strategy 
greatly  reduces  the  complexity  of  the  parameterization. 
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